Creation of macroscopic quantum superposition states by a measurement 
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We propose a novel protocol for the creation of macroscopic quantum superposition (MQS) states 
based on a measurement of a non-monotonous function of a quantum collective variable. The main 
advantage of this protocol is that it does not require switching on and off nonlinear interactions in 
the system. We predict this protocol to allow the creation of multiatom MQS by measuring the 
number of atoms coherently outcoupled from a two-component (spinor) Bose-Einstein condensate. 
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Macroscopic quantum superposition (MQS) states em- 
body the famous paradox formulated by Schrodinger [l| 
that quantum mechanics admits the existence of a cat in 
a quantum superposition of "dead" and "alive" states. 
MQS states are not only interesting from the fundamen- 
tal point of view, but are also promising for many appli- 
cations, such as precision quantum measurements 2] or 
quantum computing l3| . Most of the methods proposed 
or used so far 0, H, la 0, Hi for MQS state realization 
have been based on nonlinear unitary dynamics. 

MQS states have not yet been realized in an atomic 
Bose-Einstein condensate (BEC), although its coherently 
coupled components (in a double trap or in a spinor con- 
densate) can be isomorphous to a Josephson junction Q 
with very promising properties for the realization of MQS 
states. Yet such realization may be impeded by diffi- 
culties inherent in the proposed methods, based on the 
nonlinear dynamics of the BEC, either isolated from the 
environment . Ull [l3 ■ or perturbed by it in a specific 
("symmetrized") way |13j . 

We propose an alternative generic method of MQS 
state creation, wherein the macroscopic quantum system 
is linear, but the non-linearity is introduced by the mea- 
surement process Namely, assume that the system 
is initially in the quantum state ^'i„(a;) characterized by 
a wide spread of the macroscopic collective variable x. If 
we measure not x itself, but a certain non-monotonous 
fmiction f{x), then the state of the system collapses to 



■^out{x) = ■^„^{x)g[f{x) - fa], 



(1) 



where /o is the measurement outcome and g{y) is a 
narrow-peaked function centered at ?/ = 0. In particu- 
lar, if the Hamiltonian of the system-detector interaction 
is hPf{x), P being the detector's momentum operator, 
theng[/(a;)-/o] = (/o| exp[-iP/(x)T„]|i), where |i) and 
l/o) are the detector states before and after the measure- 
ment, respectively, and T,„ is the duration of the mea- 
surement. If the equation f{x) — /o has multiple roots, 
well separated from each other, then Eq. ([T]) describes 
a macroscopic superposition state. For the sake of def- 
initeness, assume that f{x) is even, and x — ±xo are 



the roots of the equation f{x) = /q. Then, apart from a 
normalization factor, Eq. ([T]) reduces to 



*o„t(a;) oc '^in{xo)g[f'{xo){x - xo)] + 
-^{~xo)9[f'{~xo){x + xo)]. 



(2) 



The two peaks are well-resolved if g[xof'{xo)] ^ 3(0). 
The generalization of the derivation above to the density 
matrix formalism is straightforward but cumbersome. 

The cardinal question is how to realize the neces- 
sary measurement. As a particular, experimentally rel- 
evant example, it will be shown that such a measure- 
ment is realizable on outcoupled atoms in the regime of 
coherent population trapping fis'l in a two-component 
BEC. To this end, consider a trapped BEC of atoms in 
two sublevels, |1) and |2), of the atomic ground state. 
The conjugate collective variables are the intercompo- 
nent atom-number difference and the relative phase 0, 
which are analogous to the Josephson junction conjugate 
variables, just like their counterparts in a double-well 
BEC 0, [3|. We coherently couple the levels 1 and 2 
to a common third level 0, thus forming a A-scheme of 
excitation. If the intercomponent phase is well-defined 
then the population of the state |0) is proportional to 
V^Ni + ViN2 + 2ViV2Vr^\N2COS(l), Vp being the (real) 
Rabi frequency for the /3 — transition, [see inset to Fig. 
1(a)] and Np being the mean atom number in the level 
f3, P — 1, 2. If, instead, cj) is uncertain, then we can ex- 
pect that measuring the population of the level (or its 
growth rate) will yield a certain value of cos = cos 0o 
and thereby project the state of the system to the MQS 
Eq. ^ where (p stands for x. 

The Hamiltonian in the interaction representation 
reads as 



V2(blb2 



blbo)]. 



(3) 



It is convenient to parametrize the Rabi frequencies as 
Vi = Vcosa, V2 = Vsina, V = {V^ + V^f^"^. 

Consider first the coherent dynamics generated by Eq. 
assuming for simplicity that initially (at t = 0) the 



2 



BEC is in the product of Fock states of the components 
13=1,2: 



I* 



(4) 



where |vac) is the vacuum of the atomic field and is the 
creation operator for an atom in the respective internal 
state 1/3) and in the lowest-energy motional state of the 
trap. At time t the system evolves into the state 



(6&I+li&^+eo6S)'^^|vac), 



(5) 



where = cos^ a cos Vt + sin a, ^2 = 

cos a sin a (cos — 1), Co = cosasinVt, and is 
obtained from by changing a to ^ — a. We will 
see later that accurate knowledge of A^i and N2 (or, at 
least, their difference) is essential for the detection of the 
resulting MQS state. 

Here and in what follows we assume that the number 
ng of atoms outcoupled to the level is always much 
less than Ni,N2, therefore the depletion of iVi and N2 
in the course of the evolution can be neglected. Also, to 
make the expressions less cumbersome, we assume iV2 ~ 
A'^i = N. The corresponding eigenstates, defined for the 
subsystem of atoms in the levels 1 and 2 only, are denoted 
by b)i2- Upon assuming equal Rabi frequencies, Vi = V2, 
Eq. (O reduces to 



1 



N 2j 



(-l)^(2j)!v/(2A^-2j)! 



\no)o = 



2^,t^„f^oJ-!(A^-j)V"o!(2j-no)! 

cos2J-"« Wsin"" W|j)i2|no)o, (6) 
(no!)-'/'St"°|vac). (7) 



The level populated by no atoms has a finite coher- 
ence time, caused by spontaneous relaxation, if it is op- 
tically excited, or by the translational motion of atoms, 
if it is magnetically untrapped. Most importantly, its co- 
herence is limited by the rate of the measurements that 
reveal the quantum information needed to project the 
initial state (|4]) onto a MQS. Depending on the ratio of 
the characteristic time of the evolution under Hamil- 
tonian ([3]) to the coherence time, the dynamics may be 
anywhere between two limiting regimes. If the Hamilto- 
nian evolution is much faster than the coherence time, 
we approach the limit of a fully coherent regime. The 
opposite case corresponds to the regime of continuous 
observation. Our results in the limit of continuous obser- 
vation bear similarity to those of Ref. [17] , which is a sta- 
tionary analysis of the relative phase of two independent 
condensates interfereing via a beam-splitter under ideal- 
ized conditions. Related measurement-based schemes to 
create MQS states have been suggested in Refs. 

In the limit of a fully coherent regime, a single, in- 
stantaneous measurement yields the number no of atoms 
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FIG. 1: (a) Probability distribution of the ]0) state popula- 
tion in the coherent regime. A'^i = A'2 = 1000, a = 7r/4, and 
(no) — 30. The units on this plot and the subsequent plots 
are dimensionless. Inset: schematic representation of a mea- 
surement. The sequence of detected atoms can correspond 
to either the coherent regime (overlapping signals, i.e., simul- 
taneous detection of no atoms) or to continuous observation 
(time-resolved signals), (b) Regime of continuous observation 
for the same A^i, N2, and a. Dots: {coatj)} after outcoupling 
of 30 atoms as a function of the mean time interval f be- 
tween detection events for different QMC runs. Solid line: 
approximation by r = [2IFA''(1 -I- (cos0))]~^. Inset: Four 
typical histories of establishing a definite value of (cos<jf>), as 
the number no of outcoupled atoms grows. 



in the level 0. If we drop the summation over uq in 
Eq. ([5]), we obtain (in the unnormalized form) the state 
to which the system collapses upon detecting exactly no 
atoms in the state |0). Such a measurement projects the 
state (H]) onto an eigenstate of the phase-cosine opera- 
tor cos^ = {2^/ NiN2)^^ (b\b2 + 62^1), whose spectrum is 
discrete [l^: cos0 = -1 + 2j/N, j = 0, i, 1, |, . . . , N. 
Note that terms with half-integer j's are absent in Eq. 
([6]). The crux of our method is that even if </> is not 
directly measurable, cos0 is. The corresponding condi- 
tional probability distribution of the relative phase, can 
be fairly approximated (if (j) is not too close to or ±7r) 
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by 



'\no) ^ Cph exp 



[no-iVsin^ Vt{l + cos(l))f 
2nn 



(8) 

Cph being the normalization factor. As follows from Eq. 
(I6|), the state corresponding to the double-peaked prob- 
ability distribution ([5]) is a pure state, which is therefore 
a MQS. However, such a MQS is not a sum of two Gaus- 
sian wave packets (harmonic-oscillator coherent states), 
as might be naively expected, but has far more compli- 
cated form, since terms of alternating sign appear in the 
r.h.s. of Eq. Such a behaviour precludes taking the 
limit A'^ ^ oo in Eq. ([6]). However, this can be done 
for Pph{4>\nQ), if we simultaneously set t 0, keeping 
(no) = const. 

Upon tracing out atoms in the levels 1 and 2, we ob- 
tain the probability distribution i-b("-o) of the popula- 
tion of the level 0. This distribution strongly differs 
from a Poissonian form and is plotted in Fig. 1(a). For 
TT-o 3 it is excellently fit by its quasicontinuous limit, 
whose analytic expression is cumbersome and will be 
given elsewhere. The mean value and dispersion of hq 
are, respectively, (tt-q) = N sin^ Vt and (uq) — (uq)^ — 
N sin^ Vt cos^ Vt + ^N^ sin^ Vt. 

The regime of continuous observation takes place if 

^ 7, where 7 is the inverse coherence time of the level 
0. In this case the atoms in the level are detected one 
by one. The information on the value of cos (f) is revealed 
by measuring the time intervals r between consecutive 
atom detection events. 

We have constructed a quantum Monte Carlo (QMC) 
algorithm that simulates individual measurement out- 
comes. The QMC [2l| method allows one to simulate 
the emergence of an interference pattern with a priori 
unknown phase for two spatially interfering, dissipative 
BECs, each being initially in a Fock (number) state [23 |. 

In Fig. 1(b) we present our QMC simulation of the out- 
coupling of v atoms to the level 0, starting from the state 
dl]) and detecting them one by one. The output of the nu- 
merical simulations was twofold: Firstly, we obtained the 
set of time intervals Ti, i — 1, 2, . . . , v, between consecu- 
tive detection events. Secondly, we monitored the state of 
the remaining atoms after each outcoupling event. The 
average time interval f was found to be in very good 
agreement with its estimation [2WN{1 + (cos0))]~^, 
where (cos (/>) is the mean value of the cos operator in 
the state emerging after the outcoupling of the last atom 
in the measured sequence and W = V^/j. The statisti- 
cal distribution of sequences of Ti in each QMC run was 
found to be very close to Poissonian, with the probabil- 
ity density f ~^ exp(— r/f ). This implies that the first 
few measurements of r determine quite well the value of 
(cos 0) and, hence, the subsequent evolution of the sys- 
tem [see inset to Fig. 1(b)]. 

To prove the MQS nature of the final state, one has to 



I 




FIG. 2: Left panel: Distribution of the final difference N( — 
nI as a function of the initial number difference Ni — N2 (the 
darkest area of this contour plot corresponds to the highest 
probability) for iVi « iV2 « 10^. Right panel: Probability 
distribution P of AA''i2 (final-state number difference centered 
at the initial state number difference), assuming Poissonian 
distributions of the initial atom numbers around mean values 
iVi = 1000 and N2 = 1000. The plots are labeled by the total 
r.m.s. error a of the number measurement. For both panels 
the number of outcoupled atoms 1/ = 26 and (cos (p) — 0. 



observe the interference structure in the distribution of 
the final atom-number difference, the variable conjugate 
to (f>. To this end, we plot in Fig. 2 the probability dis- 
tribution P{ANf) of atom number difference in the final 
state AA^-'' = n( - 



-N^ , as a function of the initial number 
difference AA^ = A^i — A^2- The peaks in the distribution 
of AA^-'' are separated by four atom counts, thus indicat- 
ing interference effects in the MQS state. The position 
of these interference peaks depends on the initial number 
difference, as shown in Fig. 2 (left panel). 

If we start with a mixed initial state, characterized by 
independent Poissonian fluctuations of A^i, N2 around 
their mean values iVi, N2, then the interference patterns 
can be recognized by looking at the quantity AA^i2 = 
n( — nI — (Ni — N2), the final atom- number difference 
centered at the initial number difference. This interfer- 
ence structure still persists for a measurement error of 1 
atom in determining AA^i2, cr(AA^i2) = 1.0, but disap- 
pears for cr(AA^i2) = 1.7. 

Such a high sensitivity of the MQS state detection to 
the accuracy of counting the atoms seems to be a common 
feature of measurement-based MQS phase-state creation 
methods (cf. Ref . [l9| ) . Non-demolition measurements of 
atom-number differences in high-Q cavities [23| or quan- 
tum culling techniques [2^ may help satisfy this stringent 
requirement. 

There are different possibilities to realize the A-scheme 
excitation. Preferably, the level is to be excited via 
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coupling to the ground-state sublevels by laser-induced 
single-photon transitions. In this case, to avoid compli- 
cations resulting from co-operative, bosonic-enhanced re- 
laxation [2^, we assume that Ni + N2 ^ {k^RY, where 
fci is the laser radiation wavenumber and R is the BEC 
size. The BEC should be coUisionally thin, i.e., the num- 
ber of collisions per atom moving with the momentum 
~ fikL imparted by a scattered photon should be negli- 
gible. 

Our method can be efficiently implemented with *^Rb 
condensates. The triplet and singlet s-wave scattering 
lengths are very close to each other for this isotope, and 
therefore inter- and intracomponent scattering length dif- 
ference for ®^Rb is of about 1 % [2^. This provides 
several advantages. Firstly, quantum diffusion of (j) due 
to slight difference in the mean-field interaction energies 
for the different states will destroy the intercomponent 
phase coherence on a very long time scale ~ 1 s, if we 
assume the following parameters: total number of atoms 
is about 2000, the BEC size i? « 10 /xm. Secondly, the 
inelastic losses in ^^Rb are suppressed due to the close- 
ness of the triplet and singlet scattering lengths [2^ and 
can be totally neglected for the given BEC parameters 
on time scales up to 0.1 s. Additionally, in such a small 
BEC, bosonic enhancement i2^ does not significantly af- 
fect the relaxation of the optically excited state, and the 
probability of a collision for an atom that acquires recoil 
velocity after scattering a photon is of the order of 0.01. 
Hence, a few dozens of atoms can be outcoupled without 
destroying the remaining BEC by collisions. 

The collisions with thermal atoms and scattering of 
resonance stray photons [13] destroy the macroscopic co- 
herence. In general, the rate for the MQS destruction is 
(A'^i -I- A^2) faster than the rate of excitation of an indi- 
vidual atom from the condensate (we assume that every 
such scattering is energetic enough to change the trans- 
lational state of atomic motion). Since the lifetimes of 
a BEC of the order of several tens of seconds are ex- 
perimentally feasible, an MQS can persist in a system 
with Ni K, N2 ^ 10'^ for tens of milliseconds, which is 
long enough to perform the necessary optical manipula- 
tion and detection steps whose duration lies in the sub- 
millisecond range. 

To conclude, we have proposed a novel method for cre- 
ating macroscopic superposition states in two-component 
BECs by measuring the cosine of the intercomponent 
phase. Either "snapshot" detection of the number of out- 
coupled atoms or their continuous observation have been 
shown to yield the desired result (Fig. 1). This method 
allows for fluctuations of the atom-number difference in 
the initial state, if the error in atom counting by a detec- 
tor is maintained at the level of ±1 atom (Fig. 2). Each 
measurement of cos0 yields a MQS, provided that the 
measured (j) is not too close to or ±7r. In particular, 
for the parameters of Fig. 2 (right panel), the peaks in 
two-peaked phase distribution appear to be well-resolved 



if I (cos 0)1 < 0.3. This distinguishes our method from 
the method of optical MQS generation recently devel- 
oped [i^l for homodyne measurements of quantum field 
quadratures that are suitable for photons, but not for 
atoms. A major advantage of our method compared to 
methods based on nonlinear interactions in the system 
itself [l^, [11 , l^l is that it allows one to create MQS 



rapidly enough, regardless of the nonlinearity smallness, 
well before the interaction with the environment brings 
about the collapse of the macroscopic superposition. 
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